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1. If jci, X2, • • •, ^v aie points in the unit ^-dimensional cube G.s, we may approximate to an 
integral over Gs by the average of the values of the integrand at these points. Writing 



J J ^ iT] 



(1) 



(where x\^\ x\^^^ • • • , ^i"*^ are the 5 coordinates of the point x-X we call R the ''error" of the quadrature 
formula (1). R depends on tlie integrand/and on the set of points :ri, . . . , x\: iOi R to be relatively 
small for some wide class of functions the set of points should be, in some sense, well distributed 
in the cube Gs- An appropriate sense of ''well distributed" is that, for a region A in G.s, the number 
of points of the set which lie in A should be approximately A^ times the volume of ^. Restricting 
attention to regions which are intervals with one vertex at the origin — i.e., which are of the form 

{x = {x^'\ . . . , x^''^)\0 ^ x^'^ < x\i\ i=l, 2, . . . ,5} 

for some given xo = {x{p, . . . , x^^^)eGs — 'we shall let Sn{x^o\ . . . , ^Sf^) denote the number of the 
points Xi, . . . Xn which lie in such an interval, and consider the quantity |Sa{:x|/\ . . . , x\^^) 

Let 

D^, = DUxu . . . ,:c^) = sup \S,{x^o'K . . . ,4^0-A4^ . . . , 4^')|. (2) 

xoeGs 

D% is very nearly the "discrepancy" (see, e.g., [1])^ of the set of N points. E. Hlawka [2] has 
shown that if the integrand /in (1) is a function of bounded variation and V{f) is its total variation, 
we have the bound 

\R\^V(f)^- (3) 



' Figures in bracl^ets indicate tlie literature references at tiie end of this paper. 
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Since it is very likely true that all continuous and bounded functions that one ever attempts to 
integrate numerically are of bounded variation, it is a matter of considerable interest to find sets 
of points for which Z)^ is as small as possible. 

It is easy to see that D$ ^ ^ . For if \Sn{x[^\ 1, . . . , 1) —Nx[^^ = a < 2 , then for any e > 
|Siv(:t;i^^ + €, 1, . . . , l) — N{x[^^-\-e)\ > l — a — Ne and so is > 2 if e is sufficiently small. In one 

dimension, the sequence of points Xi= ^ , i = l, 2, . . . , N, has D^= ^ for every TV. In two 

or more dimensions the situation is entirely different; K. F. Roth [3] has shown that for S > 1 
dimensions there is a positive number C{S) such that 

D* > C(S)(log N)^ (4) 



for any set of N points. It is not known that this result is best possible; there is some reason to 
guess that the exponent — ^r— may be replaceable by S— 1. In his proof Roth introduced another 
measure of unevenness of distribution, which we shall denote "En\ It is defined by 



En = En(xu . . . , xa.) = (j . . . r \Snix['^ . . . , x^^^)-Nx^'\ . . . , x^^'>\^dx^^^ . . . dx^^^'^^'' (5) 



This L^ norm is sometimes easier to handle than the maximum norm D^, and Roth in fact proved 
(4) with En in place of D^. 

In 1935 van der Corput [4] raised the question of whether there exists an infinite sequence 
^1, :x:2, . . .of points in the inverval [0, 1] such that /)*(%i, X2, . . . , ^a^) is bounded in A^. This was 
answered in the negative by Mrs. van Aardenne-Ehrenfest [5] in 1945. Roth, in [3], showed that 
this problem is equivalent to the problem of estimating a lower bound for Z)* for fixed A^ in two 
dimensions. More generally. Roth showed the existence of a constant C such that: For any A^ 
and S, if there is a set of points Xi, . . . , xn in Gs such that D^{xi, . . . , xn) = A^ there is a sequence 
x[, . . . , :r;^ in G^-i such that max D*(xJ, . . . , x^'J < C^; and conversely. It follows that for any 

n ^N 

infinite sequence x\, X2, . . • in [0, 1], /)*(xi, . . . ,xn)^ ^'(log NY^'^ for infinitely many values of A^ 
(where C is some positive constant). Van der Corput found a sequence for which he showed that 
D^ ^ log2 A^+1 (here and below "log2" denotes the logarithm to the base two); it is constructed 
as follows: for each integer ^ ^ 0, write n as 2^i -\- 2^^ -\- . . .-\-2^k with h\<h2< . . . < hk\ then 
Xn=2 (2"^^! + 2"^2+ . . . -\-2'^k). (We may picture this process as writing n as a numeral to the 
base 2 and reflecting this numeral in the "decimal point." Thus, in binary notation, %o=0, ;ci=0.1, 

1 



Xz = 0.01 , jcs == 0. 1 1 , X4 = 0.001 , etc.) Roth pointed out that then for each A^, the sequence (xu j.t 

' / 2\ ^ ^^ 

1x2, t;), . . . , i^N, 1) in G2 has the property that D^ ^ log2 A^ + 2. We shall refer to the van der 

Corput sequence as "^" and to the Roth sequence, for any A^, as "5^"'. 

J. M. Hammersley [6] and J. H. Halton [7], interested in the application of evenly distributed 
sequences to high-dimensional numerical integration, suggested higher-dimensional analogs of 
these sequences. Defining (pk{n), for any integers n ^ and /c > 1, as the number produced by 
reflecting the representation of n to the base k in the "decimal point" (more formally, if n= a\k^^ 

+ a2k^''^. . . aik^i where 0^ai<k for i=\, . . ., /, then ^k{n)=-r{axk-^'^ + . . . + aik-h)), 

K 

Halton suggested the 5-dimensional sequence defined by 

Xn = {ip2{n), (psin), (p5{n), . . ., (pp^{n)), ri = 0, 1, 2, . . . (5) 
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(where Ps is the Sth prime) and Hammersley suggested 

Xn = lj;7, (f'zin), (fsin), . . ., (^/>^, j(n)j, 7i = l,2, . . ., A^. (7) 

Halton proved that for the sequence (6), Z) * and E,w are ^ C(log A^^ and for (7) they are ^ C (log 
A0^~\ for some constants C and C. He further conjectured that the exponents K and K — l could 
be replaced by KI2 and {K— l)/2 respectively. 

2. In this paper, we shall analyze the distribution of the sequences 5^ and y in detail. Our 
results are as follows: 

Theorem 1. For the sequence 9", 

/>.?^|log2(/V) + 0(l) (8) 

and the constant 1/3 is best possible — that is, with any smaller constant (8) would not always be true. 
Theorem 2. For the sequence 9, 

^.v^^log2/V + 0(l) (9) 

and the constant 116 is best possible. 

Theorem 3. For the sequence y , 

D^ = |log,A/ + 0(l). (10) 

Theorem 4. For the sequence y\ 

£A-=|log2/V + ()(l). (11) 

From these theorems it follows that Hahon's conjecture is incorrect for the lowest-dimensional 
cases of his and Hammersley's sequences; it is very likely also to be incorrect in all dimensions. 

3. Proofs of the Theorems: 

We shall first make some notational conventions. 

If a letter is used to represent a nonnegative integer, the same letter may be used to represent 
the finite string of zeros and ones which is the binary representation of that integer. 

If q: is a finite string of zeros and ones, we shall enumerate its digits from right to left, calhng 
the rightmost one the ''zeroth digit," and write '''a = au . . . oi\ ofo." The number of digits in a 
will be denoted "||a| |"; thus ||a;c au-\ . . . ao| | = /r-h 1. For ^ i ^ k, ''a^*^" will denote the string 
ax ai-x . . . a\ ao. 

For a string a = q:a . . . ao, ".a" will denote the number -r^H — ~ h. . . + t^:^ (if a is null, 

.a will be taken to be zero). If a and /3 are strings, '"''afi''' will denote the string consisting of a and 

(3 in the order indicated; so, for example, "O.OlOo;" will denote t+q (•^)' ^^^ i^ q:=1101, this will 

4 o 

be the binary number 0.0101101, or 45/128. 

If A is any real number, " [A]'' will denote the greatest integer not greater than A, and "{A}'' 
will denote A— [A]. The distance from A to the nearest integer, which equals the lesser of {A} 
and I — {A}, will be denoted "<y^>." 

Let us now fix an integer TV > 1, and set M = [log2 N]. For any Jt;e[0,l] , let .ao«ia2 ... be 
the nonterminating binary representation of jc. Let .4 =^(x) be the string aMttM-i • • • aiOo. We 
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will first consider the sequence S^ = xo, Xi, . . ., with Xn = (p2{n). To estimate Sn{x) we proceed 
as follows: 

For Xn to be less than x, exactly one of the following M + 2 conditions must hold: 

(0'.) ao>no 

(!'.) ao=no; ai> ni 



(M'.) ao = no, ai = ni, . . ., aM-i = riM-u aM > riM 
(M+1'.) ao = no, . . ., aM-i = nM-i, aM = nM. 

These are equivalent, respectively, to the following conditions: 
(0.) Ai^0(mod2), ao=l 
(1.) n = ao(mod 2^), ai = l 
(2.) ri^ao + 2ai(mod2=^), a2 = l 



(M.) ^^ao + 2ai + . . . + 2^-iaM-i(mod 2^+i), aM= 1 
(M+1.) Ai = ao + 2ai-f. . . + 2^aM(mod 2^*+i). 



The number of integers ^ n < N satisfying condition (0.) is 
the number satisfying condition (1.) is 
the number satisfying condition (2.) is 



a2 



N-ao-2ai 



23 



etc. At most one n satisfies condition M+ 1. Thus 

"yV-(ao + 2ai + . . .-h2^"-iai-i)' 



i=0 

differs from Sn{x) by at most 1. 



2'"^ 



+ 1 



130 



Nov 



yv^=yv(|+|+ 






By the definition of M the second sum is between zero and 1. If we now set 



d=d{x, N)-- 



i=0 ^ L 



yV-(ao+. . . + 2'-V//-i) 



2'- 



-f 1- 






(12) 



then d differs from Sn{x) — Nx by at most 2; and so in proving our Theorems, for which a difference 
which is 0(1) has no effect, we may deal with d instead oi Sn{x) — Nx throughout. 

Now 

a(,-f2^/, + . . . + 2*-ia,_i. 



2^^ 



.Oai-iai-2 . . . ao 



(in binary notation), and may be written ''.OA^^ ^^";,it is equal to - (^2([2'x|). \ ;r— ^ [ may also be written 



'\mK'' Since 

r/V-(ao + . . . + 2'-V//-,) 



2' 



'2' + ^.' 



N 



we may write 



^/— V Oibi, where /;, 



+ (.A^<')-.(y<'"^)) 



l-./V^'Mf .A^<'')^.0/^('-i) 



A^ 



2/+1 



+ ]-./V<'>, 



We will first obtain a lower bound for d: hi < only if Ni = 0; therefore 



(13) 



Oibi ; 



bi 



.A"' 



(14) 



To estimate this last sum, we first note that if the string A^ (read from right to left) starts with 
a block of consecutive zeros, the .A^^'^ corresponding to those digits are zero; and so that block of 
zeros may be removed from A^. Similarly a block of ones occurring at the leftmost end of A^ may 
be removed. What remains of A^ can then be broken into disjoint substrings of the forms 01, 001, 
0001, . . . and 1, 11, 111,.... A substring 01 contributes one term to the sum — .A^^'^ for that 
/' which is the index of the zero in the substring. Since that .A^^'^ is less than V2, the contribution 
of such a substring to the sum is less than V4 the length of the substring. Similarly each block- of 
the form 001, 0001, . . . contributes to the sum a set of terms whose total is less than V4 the length 
of the block. Substrings of the form 1, 11, 111, . . . contribute nothing to the sum. It follows 
that 



d> 



M-hl 



(15) 



To obtain a sharp upper bound on d we shall require some lemmas: 

Lemma 1. For a fixed integer M ^ 1 and real number xe[0,l] , and A = A(x) as defined above, 
the maximum 0/ d(x, N) over all N between and 2^+^ — 1, is attained at N = 2A — 2^^iaM. (i.e., 
No = 0, Ni = ai_i/ori=l, 2, . . ., M.) 
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Proof: d{x, N) = ^ ai{x)bi{x, N). Changing any digit A^,- of N changes bj only for j ^ i. If a 

i=0 

change in some Ni does not change the sign of some bj,j ^ i, then it changes that bj by just ^■■^^ - 

For any i between and M, let 70 be the I'st 7 ^ i such that aj = l; it then follows that if changing Ni 
does not change the sign of bj for any 7 ^ 1, the change in aj^bj^ dominates the sum of all the changes 
in d resulting from the change in Ni. On the other hand, no bi with j ^ i changes sign unless bj^^ 
does; so if the change in Ni does change the sign of some bj, the change in aj^bj^ is still dominating. 

K 

It follows that the A^ that maximizes d must also maximize V a/6/ for every K between zero and M. 

n 

Now let rtn, ar2 - - • (n < ^2 < . . .), be the I's in ^4. To maximize V atbi we must clearly set 

i-O 

No = N\ = . . .=Nri = 0. Proceeding by induction, we assume that Ni = ai-i for i=l, 2, . . ., r^, 
and try to define Nff^, +1, • • •, ^r;^+i so as to maximize V aibi. If rk+i = rA:+ 1, we must set /Vr^.+i 

i=0 

= 1, for otherwise brj^+i would be negative. If rk+i > r/c+ 1, it is clear that ^r^^^j is maximized by 
setting A^r;^+i = l, Nrf^+2=Nrf^+-A = . . . = A^r;^.^j =0. In either case Ni = ai-i for i = rk-\-l, . . ., 
r/c+i, proving the lemma. 

Let us now, for any finite string a of zeros and ones, let 8{a) be the maximum of d{x, IV) for 
=^ A^ < 211^11+^ and A{x) = a. It remains to determine the maximum of 8(a) over all strings o: of a 
given length. 

For any integer L ^ 1, let <tl and ai be the strings lOIOl ... 01 and 0101 . . . 01011, re- 
spectively, each of length 2L -h 1 (cri = 101, (t[ = Oil). 

Lemma 2. // ^ t ^ 1, the maximum of 8{a) — t(.Q:) over all strings^ a of length 2L + 1 occurs 
at a = CTij if 1 ^ i ^ 2 the maximum occurs at a = cri. 

2 7 2 2 1 11 
Proof: We first note that 8{(Tl) = -L-\---\--——r; and that •crL = - — -—-r and .o-/'=- + 7— rr* 

3 9 9-4'^ 3 6-4^ 3 6-4^ 

The proof proceeds by induction on L. For L=l the lemma is shown true by direct calculation. 
Assuming it true for L, if a; is a string of length 2(L+ 1)+1, then q:=00^ or 01/3 or 10^ or 11/3, 
where ||/3|| = 2L+ 1. In these 4 cases respectively we have: 

3(a)-t{.a)^8{P)—^(.p) 

5(/3)-(l+^)(.i8)+ 1-^/4 

S(/3)-(| + ^)(./3)+ 1-^/2 

6(/^)-(M)(./3)-f|-f; 

we shall call these quantities 61, 62, 63, and 84. By the induction hypothesis, for ^ ^ ^ 2, 61 and 
83 are maximized by taking l3 = crL, and 82 and 84 are maximized by (3 = cri. Direct calculation 
then shows that 83 is the largest of the four when ^ i ^ 1, and 82 is largest when 1 ^ ^ ^ 2, which 
proves the lemma. 

Thus in particular the maximum of 8{a) over all strings a of length 2L + 1 is 8(cr/J, which is 

equalto||H| + 0(l). 

Furthermore, if a is a string of even length then a = 0/3 or 1^, where ||/3|| = 2L+1 for some 
integer L. In the first case 8(a) = 8{f3), and in the second 8{a) = 1 — ./3 + 8(^). Since ./3 < 1, a = 1)8 
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yields the larger 6, and then, by Lemma 2., 3{a) is maximal when 13 = ai or cri; and once again the 
maximum of 8{a) is - ||a|| + 0(l). 

Lemmas 1 and 2 together imply that for any positive integer M the maximum of d(x, N) for 
^ X ^ 1 and ^ yV< 2^+1 is ^^^^^+0(1), and Theorem L follows. 

To prove Theorem 2., we shall first show: 

Lemma 3. ( [' (Sn(x) - Nx)^ dx) "'= ^ J /|\ + 0(1). 

Proof: First of all, it is sufficient to prove the equation with d{x, N) in place oi SNix) — Nx. 
For fixed A^, audi = 0,1, . . ., M lei fiix) = ai(x)bi{x, N). Then r/U, AO =/oW +/i W + . . . +/m(^), 



an 



d 



To evaluate llie first term on the riglit of (16), we first define C,, for each /', l)y 

C,W= ^,, , . (17) 

I — .l\^'' otlierwise 

and note that Ci=hi unless ,0/^<'-') ^ ,/V<') < .O/f ■-') + ^, when C/=/>/-l. Now .0/1^'-'^ = 

-(p^iV^'x]) and so is constant on each interval (-y, . L/-= 0. 1 2'— 1: and its values on dif- 
ferent intervals differ by at least 2~'. Thus d and hj are equal everywhere except perhaps on one 
interval of length 2"', so that if we set Fi{x) = ai(x)Ci{x), we have 



|f,u)-y;w|Jxs:- and 2 '-''-l f> 

Jo ^ ,=() J ,=() J u 



2- (IB) 



Since a/ is zero on one half of eacli interval I —., —. — J and one on the oth(M- half, while C/ is constant 

r 1 1 r 1 .1 

on the interval, I Fj{x)dx = - Ci{x)dx. The values .0A^'~^^ ^ ^rr— iakrs on these 2' intervals are 

Jo 2 Jo z' 

1 2 2' 1 f 

just the numbers Wi^^ ^vTi^ ' • • ^ ^TTTT* ^^ ^^ -^^'^^T.^ C/=l-.A/<'> throughout (0,1) and F ,= 

Z Z Z Z Jo 

- (1 — .A^^'^) = - \^7—\ If -A^^'^ < -^^ vve may write .A^^'^ = /</2" ', where A is an integer between 10 and 

2' — 1. It's then easy to see that .A^^'^ < .0A^^~^^ + ^r—^ on just 2' — K of the above intervals. Then 

Ci{x)= .N^'^ on a set of measure l — KjZ' and 1 — .A"^'^ on a set of measure K/2'; so that Fi = —z- 

l / N\ Jo ^ 

= 2 \2^/ ^^^^ agam. Therefore, by (18), 

Lemma 3. will now follow if we show that the second and third terms on the right in (16) are 
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each O(log AO- For the second term this is obvious, since [/i(x)|^ 1. To bound the 3d term we 
proceed as follows: 

For ^ / ^ M, set bi{x) = - .N^^^^ aix), and yi{x)= .OA^'-^K 



Then 



n if .A^(^) ^ Yi 

[o if .A^^' < Yi 



and 



jo jo jo ^ \jo jo jo 



(20) 



(21) 



Now for any i and j, i<j, set Zi,j{x) = Yj{x) — 2^-^YiM' Recalling that yi = V2<^2([2^:e]) and 
noting that (p2{[2^x])= {2^~^(p2{[2h])} , we see that Zj, i takes on, in each interval I —., . I, the 2^~^ 
I 



2J- 



1; and each value is taken on an interval of length 2 K There- 



values 2j-[^. / = 0,1, . 

fore Zj,i, and Yi are independent functions of %, and if we define ej^{x) = efli{x) by 

n if .A^(^) ^ Zj 



1 if .A^^^) < zi 



we can conclude that 



r-*-/o4o^*=«- 



(22) 



(23) 



By reasoning as we did above about bi and C? we can see that ei{x) and ef{x) are equal except 
on a set of measure at most 2^~J'^\ and there \ef{x) — ej{x)\= 1. Therefore, by (23) 



eiej-l €i\ ej\=\\ e,(ej-e/)- eJ (e— e/)^2 
jo jo jo I IJo jo jo jo 



2J- 



From (21), therefore, 



2 ([fifj-l'fil'f)j^ 2 ^<^' 

i7lscM\Jo jo jo / ns=,vt=i/^ 



0^i<js=M 



0^i<j^M 



which completes the proof of the lemma. 

To prove Theorem 2 we must now show that 






(24) 



M / J\l \ M 

and that the constant 1/3 in (24) is best possible. We have seen that 2 (9^+1/ ^2 "^^^ {-N^^\ 

1 — .A^^^^); and we shall show that the last sum is maximized (over all ^ A^ < 2^+^) when A^ is of 
the form . . . 010101: 

For M -- 0, this is obvious. If M = A:, then A^^ OA^' or lA^', where A^' is a string with M = K-\. 

\ K-\ \ K-\ 

In the first case the sum is - (.A^')+ ^ min {M'\ \ — M^^), and in the second it is - (1 - .A^') + ^ 
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min (.N^^\ 1 — .A^^^^). The first is larger if the leftmost digit of A^' is 1, while the second is larger 
if the leftmost digit of A^' is zero; which completes the induction. For A^ of the form specified, 

X \Wi^/ ^^ equal to 1/3 log2 /V-h()(l), so (24) is established, together with the fact that the constant 

i=o \^' / 

1/3 cannot be improved upon. 

To derive Theorems 3 and 4 from Theorems 1 and 2 we will use a device of K. F. Roth [3]. 
For any A^, and ^ jc ^ 1 and ^ y ^ 1, let ai = n(y) = [Ny\. Then Sn{x, y) for the sequence^^' 
is equal to Sn{x) for the sequence 5^. Furthermore \N xy— nx\ ^ 1, so that 

sup \S'N{x,y) — Nxy\ differs from sup \Sn{x) — nx\ by no more 

< n < /V 

than 1. From Theorem 1, 

\Snix) - nx\ ^ |log2^ + 0(1) ^ |log2A^ + 0(l); 

and from the proof of Theorem 1, if ^ is the largest integer between and N whose binary rep- 
resentation is of the form 10101 ... 01, 

sup \Sn{x)-nx\ = j \\n\\ +0(1) = | ||A^|| +0(1). 

Theorem 3 follows. 

To prove Theorem 4. we use the device of Roth to replace 

(j'^J\sdx,y)-Nxyfdydyy^ by (j^(y,iy){x)-xniy)Ydxdy^ '" 
the latter, in turn, differs by only 0(1) from 



/, 



"iimf"r- 



Setting ^,(y) = (^r—j-y, we can write 

/:(i(^))*Hi/:«)-^i(/>-'^-a>f):isi/:--/:»/:«)- - 

As in the analysis of (16), we wish to show that the second and third terms on the right of (26) are 
each 0(M). For the second this is obvious; for the third we proceed so: For each i and 7 > i, set 

^j,/(k) = (]7tt[ — 2''"M— .[ \; gjj is then a function of only the M— 7th through M — i+lth digits 

of the binary representation of y, and so is independent of gi^ which is a function of the M — ith 

fi 2 

through Mth digits. Furthermore \gjj{x) — gj{x)\dx ^-^rr—; that the third term is 0{M) follows 

Jo ^ 

as in the proof of Lemma 3. 

To evaluate the first term on the right of (26), we note that 

211-570 0-66— 2 



Since (x) = { {%} ), the sum in (27) can be written as 

2^ + 1-1 



^here 



N 



2i4 



r 

2iTr 



N 



+ Eu 



0^E,= ^ 



2^ + 1 



4" 



■m 



(28) 



(29) 



Therefore 



^J„^' = ~4rl2^"l2^J + 2^^': 



(30) 



N . 2*+i 



ei. 



where |^/|< 1. 
So 



tTn J 



4 yvA 



where | fl| < 1. Since ^ 2'-> < 2''^ « A^, we can conclude that 

M 






- + 0(1) 



(31) 



(32) 



and Theorem 4. follows. 



4. An attempt to analyze the higher-dimensional sequences of Hammersley and Halton along 
the above hnes was unsuccessful, due to the greater complication of those sequences. Analogs, 
for those sequences, of Theorems 1 and 3 would be very useful in application of the sequences 
to multiple numerical quadrature. 

Intuitively, van der Corput's sequence 6^ seems to be about as evenly distributed as possible. 
This feehng, together with Theorem 1, leads us to the conjecture in section 1 that K. F. Roth's 
lower bound on /)* can be improved. 
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